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fr^ ■ tion, that proves to be at a level oia{e~^e~ — > bbcc)/a{e^e~ —>■ bb) ~ 10~^ 

for ^/s = Mz is calculated within the frames of the QCD perturbation 
theory. The cross sections for the associated production of IS'- and 2S- 

Q>^ ■ wave states of i?c-ineson in the reaction e~^e~ -^ B^bc were calculated 

fH , in the nonrelativistic model of a heavy quarkonium. The fragmentation 

f"^ • (*) 

Qh, function of 5 —> i?c is analysed in the scaling limit. The number of 

{2L|| Aftc-hyperons to be expected at LEP is estimated on the basis of the 

^ ■ assumption on quark-hadron duality. 
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Introduction. 

The study of heavy quark production processes is of great interest because this 
can give an additional check of the apphcabihty of the QCD theory to the 
description of both free and bound states. The properties of 2- and 3-jet events 
have been studied fairly thoroughly at the PETRA, KEK and LEP colliders. In 
the light of a high luminosity expected at the LEP collider, ~ 10^Z°, the study of 
other rarer processes is also interesting. Among such processes is the production 
of four quarks, bbcc, in e"'"e~annihilation. We are treating this process, on the 
one hand, as verification of the QCD theory in higher orders of perturbation 
theory, O(a^), and, on the other hand, as a further step in our understanding 
of, for example, bound states be {be) not yet discovered experimentally. The 
bound states of b- and c-quarks - i?c-mesons occupy in the mass spectrum an 
intermediate state between the families of J/ip- and T-mesons. In this respect 
the predictions of the nonrelativistic quark model describing the J/ip- and T- 
families very well can also be applicable for i?c-mesons. Therefore the discovery 
and study of S^-mesons will become a good test of QCD and QCD-inspired 
nonrelativistic potential model. The processes of i?c-meson production in e^e~- 
, pp- and also neutrino-hadron collisions have been studied in papers im - @]. 
As to the _Bc-meson spectroscopy, it was investigated in [[7|. Various aspects 
related to i?c-meson decays were considered in [§ - [p]0| . 

Using the exact formulas of QCD perturbation theory we calculate the cross- 
section a{'^e^ -^ bbec) and also the cross section for i?c-meson production in the 
process e^e~ -^ Bcbe in nonrelativistic approximation for Be- We treat these 
two processes in one work because they are closely related to each other. The 
very convenient method used to calculate the helicity amplitudes for the process 
+e~ — *> bbec can also be applied actually in the same form for calculating the 
process e~^e^ -^ Bcbe. If one considers the approximate quark-hadron duality 
to be valid then the cross section for the process """e" —* bbec can be related in 
the region of small invariant masses of a (6c)-pair with that for the production 
of bound states of i?c-mesons in the process e~^e~ — > Bcbe. Using the same as- 
sumption on quark-hadron duality one can also make certain conclusions on the 
specific value of the cross-section for the production of A^c-hyperons - a bound 
system of two heavy quarks and one light valent quark of different flavours. 

The paper is arranged as follows. The method of calculating the cross section 
for the process e~^e~ — >■ bbec is presented in Section 1. Section 2 analyzes the 
most interesting distributions for this process. The cross-section for the Bc- 
mesons production in the process e+e^ -^ Bcbe is discussed in Section 3. Section 
4 presents a brief review of the present experimental situation regarding the 
possibility of search for bound states with heavy quarks of different flavours. 



1 Calculational technique 

At the tree level, 8 Feynmann diagrams, 4 with photon and 4 with Z°-boson 
exchange in the s-channel (see Fig.l) give contribution into the process of the 
production of two pairs of heavy quarks, bb and cc: 



e+(gi) + e (^2) ^ b{pi) + 6(^2) + 0(^3) + c(p4) (1) 

The traditional techniques of calculating these diagrams by the analytical 
quadrature of the amplitude result in too cumbersome expressions. But still, 
there are some papers devoted to calculations of the cross-sections with four- 



quark final states. For example, the authors of [^ have studied within the 
frames of QCD perturbation theory the cross-section for the process gg -^ qqq'q' 
for the case of massless quarks and those of [Q have studied it also for massive 
quarks. In [^ the cross-section of the production of four quarks in e+e~- 
annihilation, e^e~ —>■ qqq'q', for the case of massless quarks has been calculated 
with the method of spinor products ||n]. As for our paper, to analyze the cross- 



section of the process e e^ -^ bbcc, where the consideration for the masses is 
obligatory, we apply a more convenient, from our viewpoint, technique of paper 
[|l^, i.e. the technique of direct numeric calculation of the amplitude. 
Let us introduce the following notations for momenta: 

gi2 = qi + q2, Pi2= - P1-P2, Psi = -ps - Pi , 

P124 = P12 - Pi , P123 = -Pl2 + P?, , P342 = P34 - P2 , P341 = -P34 + Pi , 



for currents 



Q'12 

ql2-Ml + iMzTz ' ^> 



u{pi)Yv{p; 



'2] 



P%A 

and also the following subsidiary spinors 



Ji: = ^^^^^-^^^, (4) 

P12 

Jc = -2 ' (^) 



^^124 = -^ ^JbV{Pi) , (6) 



Pl24 - ^ 



c 



{P342 +mi,) ^ 
P342 ^^6 

_. .Api2z + mc) 

M123 = U[p3)Jb^ 5- , 

_. ^ ^ (P34i + rrib) 

M34I = u{pi)Jc^ 2~ ' 

P341 - mi 



(7) 
(8) 
(9) 



In terms of the quantities introduced above the amphtude of the process 
e~^e~ -^ bbcc has the form M = J2t=i ^ii where 



Ml 
M2 
Ms 
M4 



u{Pl)[JjQ' + MVZ - aWM42 

U34i[J,Q' + Jziv'z - a'zi'Mp2) 



(10) 

(11) 

(12) 
(13) 



Here Q'^, Q'^, Q^ and w|, a|, f|, a^, f|, a^ are electromagnetic charges and 
vector and axial coupling constants of electron, h- and c-quark with Z°-boson. 

The amphtude squared should be calculated as a sum over 2® = 64 indepen- 
dent states of fermions. We choose the eigenstates of the helicity operator (Sp) 
as two independent states of the spinor ■?/': 



(Sp)V^(p, A) = \i){p,\)- 



(14) 



To calculate spinor with the specified value of 4-momentum p and helicity 
A = ±1 it is necessary to choose a specific representation. In our case we choose 
a spinor (Weyl) representation: 
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where 



1 

1 



-i 
i 



a 








a 


1 








-1 



a = {(Jx,(Jy,(yz] = 
In this representation the spinors of particles, u{p, ±), have the form 



(15) 



(16) 



u{p, +) 



( ^E+\p\ {\p\+p,) \ 

E+\p\ {Px + iPy) 



MM+Vz 



E-\p\ M+p^) 

V \JE- \v\ ij>x + iVy) ) 



:i7) 



u{p, -) 



E-\p\ i\p\+p,) 



E+\p\ {-p^ + ipy) 



V ^E+\p\{\p\+p-,) ) 
The spinors of antiparticles, f (p, ±), are defined as 

v{p,±) =Cu*{p,T) 

wliere C = i'j'^ is a charge conjugation matrix. 

The explicit form of spinors v{p, ±) is as follows: 



v{p,+] 



2\p\{\p\+p, 



( - ^E+\p\ i\p\+p,) \ 

E+\p \ {Px + iPy) 

E-\p\ {\p\+Pz) 
V \]E-\v\ {Px + iVy) ) 



v{p, -) 



2\p\{\p\+p, 



( ^ E - \P\ { -Px + iPy) \ 

E-\p\ {\p\+Pz) 

E+\p\ { -Px + iPy) 

E + \p\ m + Pz) ) 



But ii pz = — \p\ we put 



/ 



u{p, +) 







v{p, +] 



E+\p\ 

^E - \p\ 





-^E+\p\ 




u{p, -) 



v(p, -) 



V V^-N J 



V 



'^E-p \ 



-^E+\p 

y 


-^E - \p\ \ 



^/E+\p 
/ 



(19) 



(20) 



(21) 



(22) 



(23) 



The colour structure of the process e~^e~ -^ bbcc is very simple. The colour 
factor F corresponding to this process is F = (A^^ — l)/4, where N=3 is the 
number of colours. 

The amplitude of the process e+e^ -^ bbcc presented by the above formulas 
can easily be written in terms of FORTRAN codes. The resultant program 



is compact and runs fairly fast: it takes 0.03s of a VAX/VMS CPU time to 
compute the square of the amphtude summed over all polarization states. 

We have verified that our two independent programs written for calculating 
the square of the amplitude, summed over all helicity states, hold the test for 
Lorentz-invariance (boost along the beam axis) and test for azimuthal invariance 

The integration over phase space of final particles was carried out by Monte 
Carlo methods with the help of programs written for this process, that min- 
imized the spread in weights. The results of both programs coincide to an 
accuracy of the computational errors, caused by Monte Carlo methods. 

2 The cross-section for the production of bbcc 
at the Z^-peak in QCD perturbation theory 

In our numeric calculations we choose the standard set of the parameters of 
electroweak theory: Mz = 91.17 GeV, sin^ 9w = 0.23, nib = 4.7 GeV, m^ = 
1.4 GeV. The QCD coupling constant as is set equal to 

"^^^ ^ " (33-2n;)/n(QVA2)' ^^^^ 

where A = 100 MeV and for Q^ = M|,r2/ = 5 we have as{Mz) = 0.12 in 



accordance with |jT5| 



It is clear from the form of diagrams on Figs.la-d that the main contribution 
into the cross section for the process e~^e~ — > bbcc will come from the diagrams 
of Figs.lc,d because for these diagrams the virtual gluon, due to the small value 
of the c-quark mass, approaches closer the mass shell. A noticeable peak in 
the region Mcs = 2 -h Arric, related to the contribution from the diagrams of 
Figs.lc,d, is seen on Fig. 2a in the distribution over the invariant mass of cc- 
quarks, M^c = y{Pc + Pc) ■ However, the distribution over the invariant mass 

of 6fe-quarks, M^i = y{pb + Pb) ; has a completely different form (see Fig. 2b). 
A small peak in the region of small masses M^i is due to the virtual gluon 
approaching the mass shell for those diagrams in which the 66-pair is coupled 
to the c(c)-quark (see diagrams of Fig.la,b). But if the mass M^i of 6(6)-quarks 
is large, then the mass Mcs of c(c)-quarks is, as a rule, small and only diagrams 
of Figs.lc,d "work", leading to a pronounced peak in the region of large masses 
Mbi. 

The process e"'"e" -^ bbcc takes place in the second order of the QCD cou- 
pling constant asiQ"^)- Therefore it is important to make the right choice of 
the specific square of the momentum transferred Q^ in the argument as(Q^). 
As it has been mentioned above, the characteristic scale of the running q;s(<5^) 
constant may be defined by the typical gluon virtuality (M^g ~ 4m^), however, 
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Table 1: The cross-section for the process e+e~ —>■ bbcc versus the total beam 
energy ^/s for the two values of the square of the momentum transferred: Q^ = s 
and Q^ = Am^. Bracketed is the error (one standard deviation) for the last 
digital, caused by Monte Carlo method. 



Vi, GeV 


30 


50 


70 


91.17 


200 


500 


a{as{s)), pb 


0.0552(4) 


0.0782(9) 


0.192(3) 


69.7(2) 


0.061(2) 


0.014(1) 


a(as(4m^)),pb 


0.1294(9) 


0.217(3) 


0.593(9) 


234.2(7) 


0.254(8) 


0.073(5) 



the b-quark virtuality may be more greater and of the order of s, so that we 
present the results for two different choices of Q^ = {4r7ic2, s}. Table 1 presents 
the cross-sections for the process e^e'^ -^ bbcc versus the total beam energy y/s 
for these two cases. 

As is seen from this table, the difference in the definition of as leads to the 
difference in the cross-section by some times. For the most important energy 
range close to the Z'^-boson peak this difference attains a factor of 3.4. Such 
a strong dependence on a^ is related to the fact that the cross-section for the 
process e~^e~ —>■ bbcc is proportional to a^. Therefore the accurate experimental 
measurement of the cross-section for the process e~^e~ — > bbcc will allow us to 
draw a conclusion on the value of a^ in this process. 

Let us define the ratio 



Rr 



0" e ' e 



bbcc) 



ayc^e" 



bb) 



(25) 



According to our calculations, the value of Rcc is Rcc = 0.8 ■ 10^^ for Q'^ = s 
and Rcc = 2.6 ■ 10~^ for Q^ = Arn^. These values obtained from the exact 
formulas of QCD perturbation theory should be compared with the estimate 
obtained from the Monte Carlo program HERWIG (version 5.0), Rcc = (0.8 -^ 
1.8) ■ 10-2 [0]. 

Figures 3a,b present the distributions over the variables x = 2|p|/-yi, where 
\p\ is the quark momentum, for the b- and c-quark, respectively. One should pay 
attention to a completely different nature of the spectra for the b- and c-quarks. 
The 6-quarks mainly have a large momentum, whereas c-quarks a small one, 
pointing once again to the dominating nature of diagrams with fragmentation 
of 6-quarks. 

Of course, the free b- and c-quarks are not observed experimentally because 
of confinement. In fact they catch the light quark and transform into B- and 
D-mesons. To estimate the total and differential inclusive cross-section of Bc- 
meson production in pp-collisions the authors of have used phenomenological 
model of the heavy quark fragmentation into hadrons. This model describes 
quite accurately B- and D-meson production. In this approach the fragmenta- 



tion function of the i-th quark into meson, consisting of i- and j-quarks, has 



the form |17] 



Dij{x) = N x""' (1 - x)^-"^ (26) 

where aj is the leading Regge trajectory intercept, connected with the quark of 
the i-type, 7 = (1 -^ 1.5) and A^ is the normahzation coefficient: 



AT = r(2 + 7 - «, - «,)/r(i + 7 - «i)r(i - «i). 




(27) 


values of ac and ab lie in the range: ac = —(2 ^ 2.5), a^ = 


-(8- 


-9). 


jet us take q;c=-2.2, ab=-8.5 and (7 — aj) = l. Then 






Db^sS^) = 99.75 x^-' (1-x), 




(28) 


D.^B^x) = 13. U x^-^ (1-x). 




(29) 



Figures 3a,b present the spectra of B- and D-mesons, which were obtained 
by means of convolution of b- and c-quarks distributions with the fragmentation 
functions of b- and c-quarks into B- and D-mesons, respectively. 

The analysis of angular correlations also discovers a different behaviour for 
bb- and cc-pairs (see Fig.4a,b). b- and c- quarks emit mainly in opposite direc- 
tions whereas c- and c-quarks tend to be coUinear. This can also be explained 
by the aforementioned behaviour of the gluon propagator for small virtualities. 

3 The Bc-meson production cross-section 

In this Section we discuss the production of the S'-wave states of i?c-mesons. 
As to the study of the P-wave levels, we plan to make it in our next paper. 

The diagrams describing the production of i?c-niesons are obtained from 
the diagrams of Fig.la-d by combining two quark lines into mesonic one (see 
Fig.5a-d). 

The underlying assumption of our calculations is that the binding energy 
for two quarks, b, c, is much less than their masses and, hence, heavy quarks in 
the bound state B^ are actually on the mass shell. In this case the 4-momenta 
Pb and Pc of the quarks-constituents of i?c-nieson are related with 4-momentum 
P of i?c-nieson as follows: 

nib TTT-r , N 

pj, = — ^p , p^ = — £p , 30 

M M 

where M = rUb + rric is the Bc-meson mass. This assumption corresponds to 
the leading order in the effective theory of heavy quarks 



Let us also make use of the fact that the projection operators 

-7^{viP, +Mp, +) - vip, -Mp, -)} = -j={p - M)-f^ (31) 



where 



-^{viP, +Mp, +) + viP, -Mp, -)} = ^(P - M)i*{p, 0) (32) 
v{p, -)u{p, +) = -^(p - M)e*{p, +) (33) 

v{p, +)u{p, -) = -^{p - M)£*(p, -) (34) 



^^(P>0) = ^{^,p.,p„p.} (35) 

£^(p,+) = Ar{0,£,,ej„£,} (36) 

(p,-) = N{0,-s:,-e;,-el} (37) 



e^ 



and 



N = — ^ (38) 

2V2\p\{\p\ + P.) 

ex = -{\P\ + Pzf + {Px + iPyf (39) 



e 



y 



-m+Pzr + {Px+tPyy]t m 

'2{\P\ + Pz){Px + tPy) (41) 



separate the states with the specified value of the total spin S and its projection 
Sz on axis z of the system described by the spinors u{p) and v{p) {e^{p, A) is 
the polarization vector of the vector state with momentum p and helicity A). 
Then the amplitude of the process e"'"e^ -^ bbcc corresponding to the diagrams 
of Fig.5 can be expressed through these projection operators and the helicity 
amplitudes M/j^(Ai) found in Section 1 as follows: 



M(A,) = ^^ ^{0) Y: PHrnM.rHiK) , (42) 

where summation is made over the helicity states h, h of the quark and antiquark 
producing i?c-meson. The helicities of the remaining fermions are symbolically 
noted through Aj. The projection operators P^ j^ have the following explicit form 

{H=h-h)m-- 



1 

V2' 



Xk = ^l-l)'~'^'^H,o (43) 



Table 2: The cross section for the production of Be (-B*)-niesons for ^/s = Mz 
and OsiQ"^ = M^) = 0.12. Bracketed is the error (one standard deviation) in 
the last digital, caused by Monte Carlo method 



n'^+'Lj 


I'^o 


I'S, 


2' So 


2^5i 


a , pb 


0.924(1) 


1.285(1) 


0.2371(3) 


0.3168(3) 



for the "^S'n-state, 



h,h 



1 

V2' 



\H\ + -^6h,o 



(44) 



for the ^^i-state. 

Taking into account that the colour part of the i?c-nieson wave function 
is 6ij/yN the colour factor F, corresponding to the process e~^e~ -^ Bj>c, is 
F = [N^ — ly/AN'^. The value of the wave function at zero, ^1/(0), is calculated 



in nonrelativistic potential model and also in the QCD sum rules |T^ and is 
related with the decay constant fs^ of the pseudoscalar (0~) i?c-nieson and the 
constant fs* of the vector (1" 



-B*-meson as follows: 



where 



*(0) 



IBc — f. 




B* 



570 



(45) 



(46) 



The potentials of various types yield actually the same values of masses of 
lower states of i?c-niesons. For example, the mass of pseudoscalar 0~-meson 
(IS-state) is M = 6.3 GeV [0]. In this connection, when calculating the bound 
states of i?c-niesons of IS* levels we take the values of the masses of h- and 
c-quarks somewhat larger than those obtained during the production of free 
66cc-quarks and equal mb = 4.8 GeV and rric = 1.5 GeV. The mass of 25* levels 
is predicted to be M = 6.9 GeV 0. In this case the masses of b- and c-quarks 
are taken to be still larger, nib = 5.1 GeV and rric = 1.8 GeV. We would like to 
note that the authors of paper p do not take into account the effective increase 
of the masses of the b- and c- quarks, constituents of i?c-niesons, during the 
transition to higher excited states of i?c-nieson. However, this must be done 
within the frames of the "on mass shell" formalism described above. The value 
of the wave function at zero, ^(0) for 2S'-states is ^(0) = 275 MeV^^^ [0]. 

Table 2 presents the calculated cross sections for the production of Be- and 
i?*-mesons including the first excited levels. 

By analogy to what has been said in Section 2 regarding the dependence of 
the cross-section for the process e^e~ -^ bbcc on the square of the momentum 
transferred Q^ in the argument of the function ^^(Q^), similar statement for the 
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process e+e^ — *> Bj)c also leads to an essential uncertainty of the predictions 
regarding the cross-section for the i?c-niesons production. In this connection, 
the cross-sections for the process e~^e^ -^ BJjc calculated for Q"^ = s = M| 
can be interpreted as pessimistic. More optimistic predictions for the number 
of Sc-mesons to be produced can be obtained if as the square of the momen- 
tum transferred Q^ we take Q^ = 4m^. In this case the cross section for the 
production of i?c-mesons is [as(4m^)/as(s)] = 3.4 times larger. 

The above total cross section of IS"- and 2S'-states, that is equal to 2.76 
pb, is clearly a lower estimate in these calculations. Really, the total number 
of the bound states of 5c-niesons lying below the production threshold for B- 
and D-mesons is about 15. All these excited -Bc-states with probability equal 
to 1 transform into Bc-meson due to cascade radiative decays. Therefore higher 
excited states, not only IS* and 25", can make an essential contribution at least 
due to their large number. If we take into account the production of i?c-mesons, 
the above cross-section for the production of -Bc-mesons should be doubled. 

Let us define the ratio of the number of the produced i?c( -Be) -mesons, in- 
cluding also the B*{B*)-sta.tes calculated above: 

or(e+e~ -^ BJ)c) + a{e+e- -^ BM) 
a[e+e — > bo) 

According to our calculations, the value of Rb^ is Rb^ = 0.6 ■ 10^'^ for Q"^ = s 
and Rbc = 2.0 ■ 10~'^ for Q^ = 4m^. These values obtained from the formulas of 
QCD perturbation theory should also be compared with the estimate derived 
with the help of the Monte Carlo program HERWIG, Rb, = (0.1 H- 1.0) ■ 10"=^ 



le 



The authors of [^ have recently considered single production of i?c-mesons 
in e"'"e~-annihilation and obtained the analytical form for the fragmentation 
function D{x) of 6-quark into i?^*)-meson. Really, in the limit of s — *> oo, one 
can neglect the terms of the order of M2/ s anf higher. Then the differential 
cross section is equal to 

^^^ = ^"' W^^ (^'"''"(^' + ^^''"'(^) + ^3""'(^)) . (48) 

dz [1 — rz)° 

where for the pseudoscalar state one has 



9 14 

f[^\z) = 2 + {2-12r)z + {-r r'^ - 12r^) z^ + {r"^ + 2r^ + 2r^) z^ 

o o 



fPiz) = 2z\l-rz)\ (49) 

fi^\z) = 2z{l~rz) {2 + {l-6r)z + {r + 2r^)z^) , 

and for the vector state one gets 

fi^\z) = 2 + {-2-4:r)z + {3-6r + 12r^)z^ + 
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(2r - 6r^ - 4r^)z'^ + (3r^ - 2r'^ + 2r^)z^ 
fP{z) = 6z\l-rz)\ (50) 

fP{z) = 6z\l-rz) {l-2r + rz) , 

where r = rTib/Afg^Q and the calculations have been performed in the covariant 
gauge, so that the /} 2 (z) functions correspond to the squared diagrams shown 

on Figs.lc, Id, and the /3 ' (z) function corresponds to the interference of those 
diagrams. In the paper ^ the expression for the fragmentation of the 6-quark 
into the vector B* state coincides with the our one. However, the fragmentation 
function into the pseudoscalar meson is not correct in ref.0. 

Note, the scaling results (|48| - ^0|) are in a good agreement with the exact per- 
turbative calculations at y/s = rriz (see Figs. 6a, 6b). However, the nonscaling 
contributions of the order of M^/s may be essential for the both low energies 
and high mass of the meson like in the case of the 6-quark fragmentation into 
the T-particle at the Z-pole energy. 

Figures 6a,b show our fragmentation functions D(x) = ^^ oi 6-quark into 
-Bc-mesons, where x = 2\p\/y/s and |p| is the i?c- mesons momentum, for the 
production of pseudoscalar (0~) and vector (1~) i?c-niesons. The fragmentation 
function for vector mesons (see Fig. 6b) as compared with that for pseudoscalar 
ones (see Fig. 6a) has a somewhat sharper peak in the region x = 0.8 and a 
slight dip on the left of the peak. 

Using our computer program it is also easy to obtain some other distribu- 
tions. Some of them, the most interesting ones, are presented below. 

Figures 7a,b show the momentum distributions of the associative c-quark for 
the case of the production of pseudoscalar and vector _Bc-niesons. It is easy to 
understand that they are related with the distributions of Fig.6a,b by a simple 
relationship: D{xbJ — D{1 — Xc). 

Figure 8 presents the distributions in the emission angle between the pseu- 
doscalar i?c-nieson and c-quark. As is seen from this figure, the overwhelming 
number of quarks emit at small angles (< 30 — 40°) with respect to the direc- 
tion of the motion of i?c-mesons, or, as it follows from Fig.9, at small transverse 
momenta, < 5 — 10 GeV, with respect to the same direction. 

The distribution over the emission angle of 5c (-S^) -mesons with respect 
to the beam axis (see Fig. 10) can be interesting because of the effect of the 
acceptance of real detector on the number of Bc{B*)-m.esoias to be produced. 

Figs. 8-10 are shown for the case of the production of pseudoscalar (0~) 
i?c-niesons. The distributions for the case of the production of vector (1~) 
-B*-mesons have the same form, differing from the previous cases in the normal- 
ization only. 

^ In the symmetric cccc and bbbb case, when the unflavoured mesons are produced (tp, 
rjc, T, T]ii), the fragmentation function may be obtained from eqs.(Hq-pOh by the substitution 
r = 1/2. 
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4 Discussions 

So, we have calculated in our paper the cross-section for the production of 
four quarks of different flavours, bbcc, in e+e^-annihilation. At the Z^ peak it 
is 70 -7- 234 pb and depends upon the choice of the square of the transferred 
momentum in the argument asiQ"^)- When recalculated for the number of 
events, it makes up Rcc = (0.8 ~ 2.6) ■ 10^^ from the production of bb pairs or 
(1.2 -^ 3.9) ■ 10^^ from the total number of Z°-bosons. 

Since the total number of Z°-boson expected from all experiments by the 
end of 1994 is about 2 ■ 10^, this corresponds to (2.4 -^ 7.8) ■ lO'^ events with the 
production of 66cc-quarks. The experimental study of this process may be con- 
sidered as a way to verify the QCD predictions in higher orders of perturbation 
theory. 

However, the process e+e~ -^ bbcc is also interesting from another viewpoint. 
The bound states of b- and c-quarks, i.e. Be- mesons, predicted by the theory 
have not yet been discovered experimentally. Meanwhile, the discovery and 
study of i?c-mesons may yield valuable information on the behaviour of the 
quark potential in the region intermediate with respect to the families of J/ip- 
and T-mesons. Assuming the existence of quark-hadron duality, one may relate 
the cross-section for the production of the (6c)-system singlet in colour (the 
colour factor is F = (iV^ — l)^/4iV^) in the region of small invariant masses, 
Mbc, with the cross-sections for the production of the bound states of i?c-niesons 
as follows: 



Mlr.sH 



da{e^e -^ bbcc) 

be 



J dM~^ — c?^bc- = Z^^(e^e ^ B^bc) , (51) 



where mo = mh + mc< Mbc < Mb + Md + AM = Mthresh{^M ^ 0.5 ^ 1 GeV). 
Taking ttiq = 6.1 GeV as the boundary values of invariant masses and, for ex- 
ample, Mthresh = 8 GeV, we obtain the estimate for the cross-section of the 
i?c-nieson production, equal to 2.07(4) ~ 6.9(1) pb. This estimate should be 
compared with a more accurate calculation of the cross-section for the produc- 
tion of i?c-niesons obtained on the basis of the formalism expounded in Section 
3. The cross-section for the production of i?c- meson and of its first excita- 
tions, IS*- and 2S'-states, without consideration for other, (about 10), states, 
whose contribution, according to our estimates, is suppressed as compared with 
the contribution from IS*- and 2S'-states, is 2.76 -r- 9.3 pb. This is in a rea- 
sonable agreement with the value following from relationship (48) assuming 
quark-hadron duality. 

The share of events with the production of i?c-niesons is Rb^ = (0.6 -ir 2.0) ■ 
10^^ from the production of bb pairs or (0.9 -^ 3.0) ■ 10^^ from the total number 
of events. As expected, 2 ■ 10^ Z°-bosons will be produced at LEP. This means 
that the number of events with the production of iJ^-bosons, also including 
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those with the production of i?c-niesons, will be (1.8 -7- 6.0) ■ 10'^. Of course, 
the real number of experimentally reconstructed events will be less because 
one should take into account the branchings of Bc-Taeson decays into specific 
modes. As differed from J/ip- and T-mesons, i?c-mesons do not have decays of 
annihilation type into 2 or 3 gluons and decay into a weak channel. The decay 
Be — *■ J/ip + X with a branching of ~ 30% and further decay of J/ip into two 
muons with a branching of 6% is the most promising in terms of experimental 
identification of the final 30-7-110 reconstructed events with the production of 
Bc{Bc)-Taesons. 5c-mesons have quite a long lifetime and, hence, a long decay 
length, that may be helpful for identifying the vertex of Bc-meson decay. The 
availability of vertex detectors and separation of hard leptons in the jets from 
the decays of b- and c-quarks in the process e^e~ -^ Bcbc might also be helpful 
for identifying the final state, specific for the process with the production of 
Sc-mesons. However, one should take into account that this channel of B^- 
meson decay has a noticeable background coming from the decays of ordinary 
(noncharmed) rather than charmed 5-mesons into J/ip-raesons. Therefore a 
thorough analysis of this background is required for separation of the events 
with the production of i?c-mesons. 

It is interesting to note that the cross section of the production of a (bc)- 
pair with an invariant mass of a Mbc pair in the region mo < M^c < Mthresh 
(see formula (plD) is also a few ph. It is natural to assume that such a pair of 
6c-quarks, catching a light quark can fragment with a high probability into a 
colourless object, A^c-hyperon, decaying further in a cascade way. According to 
our estimates, there should be some thousands of such events per 2 ■ 10^ Z^- 
bosons. Observation of new hyperon A^c and of i?c-meson is an interesting and 
quite realistic problem of the nearer future. 
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Figure 1: Feynmann diagrams for the process e+e — > bbcc. 



Figure 2: a) Distribution over invariant mass M^c = y{Pc + Pc) of cc-pair. b) 
the same as in a) but for M^^ = \/{pb + Pi) of &6-pair. 



Figure 3: a) Distributions over x = 2\p\/^/s, where |p| is the 6-quark momentum 
(sohd hne) or B-meson momentum (dotted hne). b) The same as in a), where 
\p\ is the c-quark momentum (sohd hne) or D-meson momentum (dotted hne). 



Figure 4: a) Distribution over the angle between the directions of motions of b- 
and 6-quarks. b) The same as in a) but for c- and c-quarks. 



Figure 5: Feynmann diagrams for the process e~^e -^ Bcbc. 



Figure 6: a) Fragmentation function of pseudoscalar i?c-meson vs. the variable 
X = 2\p\/^/s, where \p\ is Sg- meson momentum, b) The same as in a) but for 
vector i?*-meson. 
. The curves are the scaling expressions (|48| - |50D 



Figure 7: a) Distribution over variable x = 2\■p\/^/s, where \p\ is c-quark mo- 
mentum for the case of the pseudoscalar i?c-meson production, b) The same as 
in a) but for the vector i?*-meson production. 



Figure 8: Distribution over the angle between the directions of motions of Bc- 
meson and c-quark. 

Figure 9: Distribution over the transverse momentum of c-quark relatively to 
the direction of motion of Sc-meson. 



Figure 10: Distribution over the angle of i?c-meson with respect to the e -beam 
direction. 



16 



